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On superconductivity in armchair carbon nanotubes. 
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We use the momentum space renormalization group to study the influence of phonons and the 
Coulomb interaction on the superconducting response function of armchair single-walled nanotubes. 
We do not find superconductivity in undoped single nanotubes. When doped, superconducting 
fluctuations can develop because of the phonons but remain small and are easily destroyed by 
the Coulomb interaction. The origin of superconductivity in ropes of nanotobes is most likely an 
intertube effect. Projections to zig-zag nanotubes indicate a more favorable disposition to super- 
conducting fluctuations. 
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Since their discovery by IjimaS in 1991, carbon 
nanotubes have attracted a lot of interest due to 
their unusual geometry and their structural and elec- 
tronic propp|ijticsa. Band calculations^ (confirmed by 
experimental'El) have stressed the one-dimensional (ID) 
character of single- wall carbon nanotubes (SWCNT). 

Recently Kociak et alu reported measurements on 
ropes of SWCNT making low-resistance contacts to non- 
superconducting (normal) metallic pads, at low voltage 
and at temperatures down to 70mK. The preparation 
technique they have used yields armchair nanotubesQ. 
Their results show signs of superconductivity below 
0.55K. The authors predict a purely electronic mech- 
anism. The question we address here is wether or not ID 
superconducting fluctuations can exist in single SWNT 
and if phonons can play any role. 

To achieve this, we perforpppperturbative renormaliza- 
tion group (RG) calculationaa'El to analyze the low energy 
behavior of a (n, n) armchair nanotube. The electronic 
and phonon parts Hq and the electron-phonon contri- 
bution He~ph of the-Hamiltonian were described in a 
previous publicatioiiE3. We added the Coulomb inter- 
action He~JtB- The electron wavefunctions come from 
a nearest-neighbor tight-binding model using a tt^ or- 
bital on each carbon atom of a sheet of graphene which 
is rolled ud in the proper way to generate an armchair 
nanotubeEll. In order to build an effective Hamiltonian 
{H — Ho + iJe_e + He^ph) for low temperatures, we dis- 
card all bands that do not intersect the Fermi level. We 
then linearize the band energies around the Fermi level, 
e~f^p{k') = (— l)'''+^pwj?fc' where k' — (k—php), 7 = 1, 2 is 
the band index, vp '^s the Fermi velocity, kp is the Fermi 
momentum, and p = ± is the sign of k. This is shown in 
Fig. 1^. The Fermi level lies exactly at the crossing point 
for a half-filled band. This is the situation we shall first 
examine. In the g-ology approach to the RG, both the 
direct electron-electron and the phonon mediated inter- 
actions give rise to twelve independent scattering ampli- 
tudes g[^J_^ and g'^^J,,^ {i = 1,2,4 and j = 1,..,4). Here 
we have adopted the notation of Krotovlla. The index i 
refers to the momentum branches (p) and j to the veloc- 
ity branches (sign of the electron velocity = {—l)'^^^p) 
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FIG. 1: Energy bands of an armchair nanotube near the 
Fermi energy. The double and single lines refer to each of the 
overlapping bands. The Fermi level is at the crossing point in 
the undoped tube and shifts when doped. 



such that 1 = interbranch backscattering, 2 = inter- 
branch forward scattering, 3 = umklapp scattering, and 4 
= intrabranch forward scattering. There are no momen- 
tum branch Umklapp processes since Akp = Sir /3a in not 
equal to a reciprocal lattice vector, where a is the tube's 
unit cell length. The coupling constants 5, were cal- 
culated by summing over all tube sites using-|a Coulomb 
interaction of the form used by Egger et aZJ13: 
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where u {u') and b {h') are related to the strength of the 
bare Coulomb interaction. The ratios {u/b,b' /h,u' /u) 
vary between (2,0.005,0.01) and (20,0.01,0.002) when 
going from a strongly screened to an unscreened inter- 
action. The phonon mediated electroivelectron coupling 
constants in the non-adiabatic regimell^, gp^^, are given 
byE: 
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where ge-ph.-yil) is the electron-phonon interaction with 
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each of the two electrons involved and ujq is the phonon 
frequency. One has q~2kp for backward scattering and 
g « for forward scattering. These parameters are valid 
for temperatures or energies {ks — 1) smaller than the 
phonon energy. _ 

The symmetry of the electronic functions liil leads to 



the following sign constraints for the g. 
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We then perform a one loop RG calculation, which is a 
generalization of the one band calculationElJ, by consider- 
ing only the diagrams presenting a logarithmic divergence 
with temperature, that is with electrons on different ve- 
locity branches. The result of this procedure yields the 
same flow equations foe, the coupling constants that were 
found by Krolpv et al£3. At half-filling, we use a two cut- 
off approaclo with W > fnoph where W is the bandwidth 
and the renormalization starting energy scale while fiLOph 
is the phonon energy. For energy (or temperature) scales 
Eq = We~^ > hujph, the phonons are adiabatic and the 
9e-ph,-y renormalize only in the random phase approx- 
imation. When Eq reaches fkOph, the phonons become 
non-adiabatic, the 9^''pi^ become active and are added to 

the gl e-e- Contrary to the usual case of a single band 
crossing the Fermi level, the coupling constants in our 
case scale towards the strong coupling sector. An an- 
alytical solution of the renormalization equations seems 
to be impossible in the general case, which forces us to 
resort to a numerical solution. 

The important physical properties of the system can 
be probed through its various response functions. These 
have the coupling constants as input and measure the 
relative importance of the underlying fluctuations. It is 
important to realize that, because of the interband inter- 
action terms ff2^\ 54''^), the particle- hole pairs 
in the Peierls channel and the particle-particle pairs in 
the Cooper channel will be evolving in both bands. It is 
thus necessary to define the following response functions 
in Matsubara-Fourier space 



Xn 



■Mtx 



with K 



±, M 



drdx e 

x(O;^^-^(x,r)to;-^'^(0,0)), (4) 
±, and q = (g,a;,„). The Fourier 



transforms OJ^' (q) are defined by 

0;-''\q) = -^(O^iq) + MO;iq)^). (5) 

In the Peierls channel, one defines 
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in which fi — Q stands for charge density (CDW) opera- 
tors and ^ — 1,2, 3, for spin density (SDW) ones. Here 
(To and cri^2,3 are the identity and the x, y, z Pauli matri- 
ces, respectively, and ip-y^p^dk) annihilates an electron of 
spin a in band 7 having momentum k in branch p. We 
have introduced the parameter M by which we differen- 
ciate between on-site and bond-order for charge and spin 
correlation functions as defined in Ref. |l^. In the Cooper 
channel, one has 
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where /i = are singlet superconducting (SS) operators 
and ^ — 1,2,3 are triplet superconducting (ST) ones. It 
is through these band-entangled operators that our re- 
sponse, functions are different from the ones of Krotov 
et alt3 who only used the untangled (0+ ± 0~ ) opera- 
tors. Our definition of the response functions leads to a 
fundamentally different behavior with temperature. 

In order to calculate the evolution of the response func- 
tions with the energy scale Eq (or temperature), we in- 
troduce the auxiliary response functions x^''^ defined 
through 



.Ml 



A) = - 



1 



(10) 



where I = ln(iy/_Eo). We deduce the following renormal- 
ization equations 
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FIG. 2: Typical flow diagram for the response functions of a 
doped armchair nanotube. 

with g,p^ = g-^^/TrvF- CDW' and SDW' refer to the 
anomalous q ^ interband situations while SS and TS 
refer to the q ~ 2kp ones. 

We first report on the calculations without the 
Coulomb interaction, that is solely with the phonon- 
mediated effective electron-electron interactions (at tem- 
peratures below the Debye temperature for non-adiabatic 
interactionsEJ). This was done for arbitrary amplitudes 
but with the sign constraints given in Eq. (|3|). We find no 
sign of a dominant superconducting response in a single- 
walled armchair carbon nanotube. Charge correlations 

iXcDW ^'^'^ XcDw) ^'^^ ™^ cases the most diver- 
gent and open a pseudo-gap which subdues the super- 
conducting fluctuations. The introduction of Coulomb 
interactions, through u and 5, reinforces this tendency 
even more. Kociak et a/.B, however, mention the possi- 
bility that the band occupancy might not be exactly 1/2. 
If this were the far as we can estimate from a 

tight-binding calculation, the sign constraints mentioned 
in Eq.( ||) still hold. However, the Fermi level would 
shift by Ai? = dvphp/^ {d is the doping level) and the 
two bands would have different Fermi momenta. As a 
consequence, for energy scales below AE the gf^ would 
vanish because of longitudinal momentum conservation. 
Moreover, the interband backward scattering g\ ^i^, which 
was previously zero (see EqJ^), now is ~ d^gi- We now 
use a three cutoff procedure W > hujph > AE, which 
applies up to 8% doping level where AE w hcuph- In 
the case where only the phonon mediated interactions 
are considered we have estimatedcJ gi — — 0.3/n and 
52 = 54 = — 0.1/n. Superconducting correlations Xs's 
are found to dominate for i (= ln(W/r))greater than a 
critical £s {— ln(Vl^/T5)). This is shown for a typical run 
in Fig. H (only those giving the largest enhancement are 
shown). Fig. |^ shows the crossover temperatures Ts be- 
low which this occurs for various doping levels and differ- 
ent diameter SWCNTs . The turning on of the Coulomb 
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FIG. 3: Temperature at which the superconducting fluctu- 
ations dominate as a function of the doping level for three 
(n, n) armchair tubes. Only phonon-mediated interactions 
are considered. 



interaction lowers Ts and eventually destroys the super- 
conducting fluctuations dominance for u > |(7i|. At this 
point, we can conclude that superconducting fluctuations 
in single armchair nanotubes originating from a phonon 
mechanism are possible at an appreciable doping level 
provided the Coulomb interaction is very small. This 
last condition seems restrictive in view of our own eval- 
uation u >> \gi\ and th.e evidence of strong electronic 
correlations in SWCNTlia. Moreover, the finite length of 
the najiotubes results in a discretization of the energy 
levelalj. Consequently, the energy scale in the renormal- 
ization group procedure cannot be smaller than this spac- 
ing. We estimate this would occur at best at IK which 
is still laijger than the temperatures at which the work 
of Kociaka shows any indication of superconductivity It 
thus seems likely that the superconductivity seen in ropes 
takes place by a coupling effect between the nanotubes. 
At the dimensionality crossover temperature Ta; ~ t^/n, 
where t± is the net inter-tube hooping amplitude, inter- 
tube hopping becomes coherentEJ. Any fluctuating su- 
perconducting pair that might exist will thereafter be 
able to coherently tunnel between tubes. This surely oc- 
curs much before the SWCNT discrete spectrum is felt. 
The single tube superconducting fluctuations need not 
be dominant for this to take place. The existence of frus- 
tration in ropes of close packed SWCNT would prevent 
the further development of bond-order wave deforma- 
tions (associated with a Kekule CDW modulatior£3) or 
magnetic order and allow superconductivity to develop. 
Moreover, the interchain particle-hole interactions that 
might exist at could also play a role in enhancing su- 
perconductivity through a mechanism similar to the one 
proposed for organic materialslilrEj. All of this is consis- 
tent with the small temperatures observed for supercon- 
ductivity in ropes of armchair nanotubes. 

The recent discovery of superconducting fluctuations 
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at 15 K in a single (5, 0) zigzag thus seems 

puzzling in view of the above analysis of armchair nan- 
otubes. In this specific instance, the very strong a — ir hy- 
bridization due to the very small diameter of the SWCNT 
changes the band structure in a dramatic way and makes 
the tubes, metallic instead of being the expected insulat- 
ing stataa. The band structure might share siruilarities 
with the one in a conducting (6, 0) zigzag SWNTE3. The 
band structure in the vicinity of the Fermi level would 
then show the crossing between a non-degenerate band 
(1) and doubly-degenerate bands (2, 3) of opposite cur- 
vature and having a different transverse angular momen- 
tum. Because of the twofold degeneracy, the Fermi level 
no longer lies at the crossing point but is offset so that 
the bands (1) and (2, 3) have different Fermi momenta. 
Moreover, the tight symmetry relationship of the arm- 
chair bands exists no longer between (1) and (2, 3). There 
are also many different phonons contributing to effective 
phonon-mediated electron-electron interactions. Finally, 

-(3) 

all gl will vanish because of longitudinal momentum 
conservation. This, we believe, is quite sufficient to allow 
for important superconducting fluctuations in the zigzag 
tubes. Purely as an illustration and only for crude or- 
der of magnitude estimates, we used the model deV|elp 
oped above with only two cutoffs W — AE ~ 0.3eVc3, 



hujph '■'^ 0.166 eV, and n = 3 to account for the smaller 
tube diameter such that gi — —0.1, 52 — —0.03. We 
find a dominance of the superconducting response below 
room temperature with just the phonons. This is much 
larger than for the armchair tubes and might explain 
the origin of the superconducting fluctuations observed 
by Tang et alM. But again adding a Coulomb interac- 
tion quickly reduces this temperature. Superconductiv- 
ity disappears again for u > \gi \ ■ The results of Tang 
would thus indicate an unexpectedly small Coulomb in- 
teraction in zig-zag nanotubes. A more detailed analysis 
of the three-band RG will be given elsewhere. 
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